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The role of system size and momentum dependent effects are analyzed in multifragmenation by 
simulating symmetric reactions of Ca+Ca, Ni+Ni, Nb+Nb, Xe+Xe, Er+Er, Au+Au, and U+U at 
incident energies between 50 MeV/nucleon and 1000 MeV/nucleon and over full impact parameter 
zones. Our detailed study reveals that there exist a system size dependence when reaction is simu- 
lated with momentum dependent interactions. This dependence exhibits a mass power law behavior. 
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I. INTRODUCTION 



The last gem in the field of heavy-ion collisions, 
namely, multifragmentation has always attracted the- 
oreticians as well as experimentalists [IHEl- Primarily 
due to the several hidden phenomena that need deeper 
investigations and secondary due to its connection with 
nuclear equation of state - a question which always 
has captured a central place in the research of nuclear 
physics. The knowledge of the nuclear compressibility is 
not only relevant for the nuclear physics, it is also vital 
for other branches such as astrophysics. One should, 
however, note that the compressibility depends not only 
on the density but also the entire momentum plane. In 
other words, equation of state apart from the population 
of nucleons also depends upon their relative velocities. 
This can also been seen from the optical potential where 
strong momentum dependence was reported in the 
literature @. 

The momentum dependence of the nuclear equation 
of state has been reported to affect the collective flow 
and particle production drastically 0, 0-Qjl- Some 
initial investigations also point toward its important 
role in multifragmentation [H, |{| [TH, [I?}- Due to the 
repulsive nature of momentum dependent interactions 
(MDI), nuclei propagating with MDI are reported to 
emit nucleons and light complex fragments. However, 
the rms radii of single nuclei is not affected significantly 
Q. One has to keep in the mind that the response of 
momentum dependent interactions also depends on the 
system size. For example, it has been shown by Sood 
and Puri [Toj that momentum dependent interactions 
push energies of vanishing flow to significant lower levels 
for 12 C + 12 C system, whereas for heavier systems, the 



trend is just opposite. In a recent study, Singh and Puri 
fllj predicted a power law for the system size effects 
in multifragmentation. In their study, a simple static 
equation of state was used. In view of the above facts, 
it is challenging to investigate the role of momentum de- 
pendent interactions with respect to multifragmentation 
& see how system size affects the outcome. Our present 
aim, therefore, is three fold atleast: 

(i) to understand the role of momentum dependent 
interactions in multiframcntation, 

(ii) to study the system size effects in the presence of 
momentum dependent interactions and (iii) to find a 
scaling to these system size effects. 

This study is done within the frame work of quantum 
molecular dynamics model. The section II deals with 
the model, section III discuss the results. Our results 
are summarized in section IV. 



II. THE MODEL 



The QMD model is a time dependent many-body the- 
oretical approach which is based on the molecular dy- 
namics picture that treats nuclear correlations explicitly. 
The two dynamical ingredients of the model are the den- 
sity dependent mean field and the in-medium nucleon- 
nucleon cross-section [l8j]. In the QMD model, each nu- 
cleon is represented by a Gaussian wave packet charac- 
terized by the time dependent parameters in space fl(t) 
and in momentum pi (t) [l[ . This wave packet can be rep- 
resented as: 
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The total n-body wave function is assumed to be a direct 
product of the form: 



(2) 



The model uses its classical analog in term of Wigner 
function [H. 



fi(r,p,t) 
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(3) 

The parameter L is related to the extension of the wave 
packet in phase space. This parameter is, however, kept 
fixed in the present study. The centroid of each nucleon 
propagates under the classical equations of motion [l|: 



dfl dH dpi dH 
dt dpi dt dfl 

The H referring to the Hamiltonian reads as: 
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(EOS) with comressibility K = 200.MW. A standard 
energy dependent nucleon -nucleon cross-section due 
to Cugnon was also used In brief, we followed the 
time evolution till the end of the reaction which, in the 
present study, is 200 fm/c. 
Here each of the reaction ^Ca +20 Ca, 

93 /V/)-|_93 Afh 131 -v> 1 131 y„ 167 

and g2 8 C^+|2 8 U was simulated for 100 events at incident 
energies between 50 and 1000 MeV/nucleon using 
different collision geometries. By using the symmetric 
(colliding) nuclei, system size effects can be analyzed 
without varying the asymmetry (and excitation energy) 
of the system. It is worth mentioning that the exper- 
imental studies by the MSU mini ball and ALADIN 
@i Hil groups, varied the asymmetry of the reaction 



whereas plastic ball [21[ and FOPI experiments [22 1 
are performed for the symmetric colliding nuclei only. 
In the following, we first discuss the time evolution of 
different reactions and then shall address the question 
of momentum dependent interactions and system size 
effects. 



Our total interaction potential V tot is a composite of var- 
ious terms: 
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Here m = 1.5 fm, t 3 = -6.66 MeV, t 4 = 1.57 MeV and 
t 5 = 5 x 10" 4 MeV" 2 . 

The momentum dependent interactions can be incorpo- 
rated by parameterizing the momentum dependence of 
the real part of optical potential @ . 



III. RESULTS AND DISCUSSION 

We simulated each reaction at various time steps and 
stored the phase space. This phase space needs to be 
clustcrized. In the present study, we clusterized the phase 
space using the minimum spanning tree (MST) method 
which binds two nucleons in a fragment if their centroids 
are closer than 4 fm. In other words, we demand: 



\n -fj\ < 4 fm 



(10) 



At present, various reactions were simulated with soft 
and soft momentum dependent (MDI) equations of state 



A. Time evolution 

The density of the environment surrounding the nu- 
cleons of a fragment plays crucial role in deciding the 
physical process behind their formation. In figlU we dis- 
play the average density < p/po > reached in a typical 
reaction as a function of the time. The average nucleonic 
density < p/po > is calculated as (l8j 



<p/p > = 
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with fl and fj being the position coordinates of i th and 
j th nucleons, respectively. We here display the average 
density at incident energies of 50 and 400 MeV/nucleon. 
In addition, two colliding geometries corresponding to 
6 = and b = 0.6 are also taken. The reaction 
(at low incident energies) preserves most of the initial 
correlations and hence only a small change in the density 
profile occurs. This trend turns to a sharp decrease 
at higher incident energies. This is due to the fact 
that higher incident energies lead to highly unstable 
compressed zone which, does not sustain for a long time 
and as a result fast emission of nucleons occurs. 
A similar trend is also seen for the collision profile 
which has a direct relation with the density reached in a 
reaction. 

Let us now turn towards multifragmcntation. In figlU 
we display the heaviest fragment A max survived in a 
reaction. The medium mass fragments (MMF's), defined 
as 5 < A < 9, are displayed in fig 01 The top panel in 
both figures is for 50 MeV/nucleon, whereas bottom 
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FIG. 1: (Color online) Average density (p/po) as a function 
of the time. The top panel is at 50 MeV/nucleon, where bot- 
tom panel represents the reaction at 400 MeV/nucleon. The 
left and right hand sides represent, respectively, the central 
collision 6 = and peripheral collision b = 0.6. All reac- 
tions represent symmetric colliding nuclei X + X, where X 
represents the reacting elements. 



FIG. 2: (Color online) Same as figH] but for the time evolu- 
tion of the heaviest fragment A max as a function of the time. 
The shaded area corresponds to density higher than normal 
nuclear matter density po- 



panel is at 400 MeV/nucleon. The time evolution of the 
fragments reveals many interesting points: The heaviest 
fragment A max survived in the reaction of heavier 
systems struggles for a longer time. The soft momen- 
tum dependent (SMD) equation of state destroy most 
of the nucleon-nucleon correlations, therefore, A max 
obtained with SMD EOS is lighter than corresponding 
soft equation of state Consequently, there is an 

enhanced emission of the free nucleons, LMF's (light 
mass fragments) and MMF's. We also see an appreciable 
enhancement in the nucleonic emission (not shown here) 
at all incident energies and impact parameters. For 
heavier systems, emission of nucleons continues till the 
end of the reaction. This is due to the finite collisions 
happening at the later stage as well as due to the longer 
reaction time at these incident energies. 
It takes longer time for A max in heavy systems to be 
stabilized compared to lighter nuclei, where saturation 
time for A max is much less. The excited A max in heavier 
system continues to emit the nucleons till the end of the 
reaction. This time is, however, much shorter in lighter 
nuclei. 

The multiplicity of MMF's has a different story to tell. 
We now see more fragments in central collisions at 50 
MeV/nucleon compared to peripheral collisions. As we 
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FIG. 3: (Color online) Same as in fig[T] but for the time 
evolution of multiplicity of MMF's. 
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FIG. 4: (Color online) Final state multiplicity of the medium 
mass fragments per projectile nucleon as a function of pro- 
jectile mass Ap. The top, middle and bottom panel at left 
side represent, respectively, the reactions at 50, 200, 600 
MeV/nucleon. The corresponding right hand side represent, 
respectively, the reactions at 100, 400, 1000 MeV/nucleon . 
All curves are fitted using function Y = CA P . 
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FIG. 5: (Color online) Parameter r (appearing in the power 
law function Ap) as a function of incident energy. The panel 
displays the value of r for medium mass fragments. 



increase the energy, the MMF's production decreases 
considerably. This is valid for the central collisions only. 
The peripheral collisions yield almost same MMF's. As 
noted in ref. [ll|, the static soft EOS is not able to break 
the initial correlations among the nucleons in peripheral 
collisions. As soon as momentum dependence is taken 
into account, the initial correlations among the nucleons 
are destroyed, resulting in large number of MMF's. The 
simple static EOS fails to transfer the energy from the 
participant to spectator matter. In other words, MDI 
suppresses the production of free nucleons and LMF's 
while the production of MMF's is enhanced. Overall, we 
observe enhancement in the multiplicity of medium mass 
fragments with MDI compared to static equation of state. 



B. Final state fragment distribution 

As we know, measurements are always done at the 
end of the reaction. The reaction time is chosen to be 
t = 200 fm/c. It is based on the fact that directed flow 
saturates by this time. Therefore, it will be of interest to 
see whether the final state fragment distribution has mo- 
mentum and system size dependence or not. We display 



in fig HI the reduced multiplicity (multiplicity /nucleon) 
of medium mass fragments (MMF's). In experimental 
observations (e.g. FOPI and ALADIN), the nuclear 
matter is divided into spectator and participant parts 
23]. In our calculations, this division is made by 
splitting the reaction into different impact parameter 
zones that can be related with the spectators/participant 
matter. The top panel in figU displays the multiplicity 
of MMF's at 50 and 100 MeV/nucleon, whereas bottom 
panel is at 600 and 1000 MeV/nucleon, respectively. 
The middle panel represents the outcomes at 200 and 
400 MeV/nucleon. Each window of the panel contains 
four different curves that correspond, respectively, to the 
scaled impact parameter values of b = 0.0,0.3,0.6,0.9. 
The wide range of incident energy between 50 and 1000 
MeV/nucleon and impact parameter between zero and 
b max gives opportunity to study the different dynamics 
emerging at low, intermediate and higher energies. 

The energy received by the target in peripheral colli- 
sions is not enough to excite the matter above the Fermi 
level, resulting in fewer light fragments. As a result, 
emission of heavier fragments takes place. Irrespective 
of the incident energy and impact parameter, the mul- 
tiplicity of medium mass fragments (5 < A < 9) scales 
with the size of the projectile that can be parametrized 
by a power law of the form C(Ap) T ; is the mass of 
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the projectile. The values of C and r depends on the 
size of fragments as well as on the incident energy and 
impact parameter of the reaction. 

In fig|5j we plot the power law parameter r as a function 
of incident energy and impact parameter. Although, no 
unique dependence occurs for t, we can correlate some 
of its values. No physical correlation can be extracted 
for the central collisions. This is perhaps due to the 
complete destruction of initial correlations, moreover, 
as a result, even no collective flow has been observed in 
central collisions [241 ]. 



1000 MeV/nucleon and over full geometrical overlap. 
The inclusion of momentum dependent interactions 
leads to less freeze out density in all colliding systems. 
This happens due to repulsive nature of momentum 
dependent interactions. The system size effects are 
found to vary with reaction parameters and incident 
energies. The multiplicity of medium mass fragments 
can be parametrized in term of a power law. This is 
true for a wide range of impact parameters and incident 
energies considered here. However, the parameter r does 
not have unique value. 



IV. CONCLUSION 

We have studied the role of momentum dependent in- 
teractions in fragmentation by systematically analyzing 
various reactions at incident energies between 50 and 
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